Hot electron driven enhancement of spin-lattice coupling in 4f ferromagnets observed 
by femtosecond x-ray magnetic circular dichroism 
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Femtosecond x-ray magnetic circular dichroism was used to study the time-dependent magnetic 
moment of 4/ electrons in the ferromagnets Gd and Tb, which are known for their different spin- 
lattice coupling. We observe a two-step demagnetization with an ultrafast demagnetization time of 
750 fs identical for both systems and slower times which differ sizeably with 40 ps for Gd and 8 ps 
for Tb. We conclude that spin- lattice coupling in the electronically excited state is enhanced up to 
orders of magnitude compared to equilibrium. 

PACS numbers: 78.47.J-, 71.38.-k, 75.70.Ak, 78.70.Dm 



Laser-induced magnetization dynamics has high po- 
tential for ultrafast data-storage applications [1] and a 
microscopic understanding of the underlying processes 
is essential for device optimization and tuning. In this 
context switching the magnetic order by intense, ultra- 
short laser pulses explores the speed limit of magnetic 
recording. Next to its technological relevance magneti- 
zation dynamics driven by femtosecond (fs) laser pulses 
challenges our microscopic understanding of magnetism: 
(i) Bigot et al. [2] and Zhang et al. [3] suggest that 
the light field is involved in magnetization dynamics, (ii) 
Battiato et al. propose superdiffusive spin transport as a 
mechanism of ultrafast demagnetization 0. (iii) Koop- 
mans and coworkers have developed an empirical model 
based on spin-orbit mediated electron spin-flip scattering. 
Their concept implies a material dependent demagnetiza- 
tion time and connects itinerant and rare earth ferromag- 
nets [5] . Ultrafast laser-induced magnetization dynamics 
has been established for the 3d metals and a number of 
alloys 

In view of angular momentum conservation 
a change in the magnetization M requires transfer of an- 
gular momentum from M to some other reservoir. The 
crystal lattice is a prominent candidate here, which turns 
spin-lattice coupling into an essential, but barely investi- 
gated interaction in ultrafast magnetization dynamics. 

In this letter we report on laser-induced magnetization 
dynamics in the lanthanide ferromagnets Gd and Tb. By 
time-resolved x-ray magnetic circular dichroism (XMCD) 
at the M5 absorption edges we probe directly the 4/ 
magnetic moment, out of reach for magneto-optical tech- 
niques. We identify for both materials two separate de- 
magnetization processes, a slower quasi-equilibrium one 
and an ultrafast one active in the electronically excited 
state. The time constants for the slower process differ for 
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FIG. 1: (color online) Orbital wave-function distributions 
within an hep unit cell for L, m — and 3. The m — 3 
non-spherical distribution of Tb couples to the ion cores via 
single ion anisotropy, which is absent for the spherical m = 
state of Gd. 



the strong direct spin-lattice coupling in Tb (8 ps) and 
the weaker indirect interaction in Gd (40 ps). The ultra- 
fast process agrees for both elements (0.74 vs. 0.76 ps) 
and is active while hot electrons are present. It in- 
volves an enhancement of the indirect spin-lattice cou- 
pling, which leads to a pronounced increase in the mo- 
mentum transfer rates from the magnetization to the lat- 
tice in Gd by as much as 50 times. 

The heavy lanthanides Gd (4/ 7 ) and Tb (4/ 8 ) are well 
known for their magnetic properties as a function of oc- 
cupation of the 4/ orbital. While the spin quantum num- 
ber S decreases as the 4/ shell is more than half filled 
(Gd S = 7/2, Tb 6/2), the orbital quantum number L 
increases (Gd L = 0, Tb 3) [15]. The magnetic mo- 
ment per atom [i at follows Hund's rules (Gd 7.55 /ig, Tb 
9.34 fiB [l6[), where the excess from the integer value is 
attributed to spin polarization of the 5d6s valence elec- 
trons. Fig. [1] depicts the L, m = and L, m = 3 angular 
distribution of the 4/ orbital of Gd and Tb, respectively; 
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FIG. 2: (color online) a) X-ray transmission at the M5 absorp- 
tion edges of Gd and Tb films recorded for opposite magne- 
tization direction (black and gray lines) with 10 ps circularly 
polarized x-ray pulses, b) XMCD signals of Gd and Tb before 
and 200 ps after laser excitation (solid and dotted lines). 



m is the magnetic quantum number. A pronounced cou- 
pling of the orientation of n at to the neighbouring ion 
cores and hence to the lattice follows for Tb from the 
non-spherical 4/ distribution since spin-orbit interaction 
couples the direction of the spin moment to the 4/ or- 
bital. Such a non-spherical distribution links a rocking 
of the atomic magnetic moment directly to a lattice vi- 
bration and vice versa. For the spherical distribution of 
the half filled Gd 4/ shell this direct coupling is absent 
(Fig. [I]). Indeed, the magnetic anisotropy constant K2 
describing the energy required to rotate M with respect 
to the basal plane of the hep lattice is in Gd more than 
two orders of magnitude smaller than in Tb [15] . Also 
magnon excitations reflect this difference in L. Avoided 
crossings in the magnon dispersion of Tb explaine d by 
magnon-phonon coupling [l7| are absent in Gd [HI, [l9j • 
The magnetic anisotropy in Gd is, however, non-zero due 
to 4/ — 5d coupling and the spin-orbit interaction of 5d 
electrons [20]. We refer to such a valence electron medi- 
ated spin-lattice coupling as indirect. 

Optical pump - x-ray probe experiments were per- 
formed at the femtosecond slicing facility of BESSY 
II 



21] . The 5d6s valence electrons were excited by 



1.5 eV laser pulses of 50 fs duration at a fluence of 
F = 3 — 5 mJ/cm 2 with the sample held in an ap- 
plied magnetic field of 5 kOe at an equilibrium tempera- 
ture of 140 K. We measured x-ray transmission for poly- 
crystalline Y(50 nm)/R(10 nm)/Y(5 nm) films grown on 
a free-standing 0.5 fim thick Al substrate; R = Gd,Tb. 
The x-ray photon energy was tuned to resonantly ex- 
cite the 3c?5/2 core- level electrons to the unoccupied 4/^ 
states with a binding energy of 4 eV above [22| . Since 
optical transitions between 4/ and 5d require photon en- 
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FIG. 3: (color online) Time-dependent XMCD signals for Gd 
(top) and Tb (bottom) measured by 10 ps x-ray probe and 
50 fs laser pump pulses. Solid lines indicate fits to the data. 
The inset depicts Gd data in a smaller time window with the 
actual time-resolution of 16 ps indicated. The biexponential 
fit (solid line) highlights the two step demagnetization process 
and the dashed line indicates the behavior expected for an 
instantaneous first step. 



ergies far above 1.5 eV, 4/ levels do not participate in 
the optical excitation [23[ and can be used as a reliable 
monitor of M [241 ]: pump-induced refilling of 4/ levels 
and saturation effects do not affect our measurements. 

Figure [2^ shows the transmission spectra for M par- 
allel (+) and antiparallel (— ) to the helicity of circularly 
polarized x-ray pulses without laser excitation. XMCD is 
determined from the difference of the absorption for op- 
posite M. Comparing XMCD signals before and 200 ps 
after laser excitation (Fig. [2)3) exhibits a pronounced 
pump-induced change. The sum of the spectra (not 
shown) remains unaffected even though the temperature 
is increased by the optical excitation. This guarantees 
that the change in XMCD is a purely magnetic effect. 

We proceed to the magnetization dynamics and ana- 
lyze the time-dependence of the XMCD signal. At first, 
we employed x-ray pulses of about 10 ps duration, which 
are available in the low-a operation mode of BESSY II 
[25] . Fig. [3] depicts the time-dependent XMCD signal for 
Gd and Tb normalized to the value before optical exci- 
tation. For both materials we find a pronounced demag- 
netization, but the detailed behavior is different. For Gd 
the minimum of M is reached after 200 ps in a two-step 
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process. The inset indicates that about half of the final 
demagnetization occurs within the 10 ps pulse duration 
of the x-ray pulse, while the second process lowers M un- 
til 200 ps. We fit the Gd data by a biexponential decay 
convoluted with the x-ray pulse duration and determine 
a characteristic time constant of r^ d = 40 ± 10 ps for 
the slower process. Electrons and phonons have equi- 
librated after 1 ps [26j. Therefore, we refer to delays 
> 1 ps as a quasi-equilibrium and < 1 ps as an elec- 
tronically excited state. The obtained r e ^ d is character- 
istic for the weak indirect spin-lattice coupling in Gd 
(L = 0, c.f. Fig. [T]) in quasi-equilibrium, since it is 
^>1 ps. Our finding substantiates previous experimen- 
tal and theoretical results [27H29j. From the change in 
M stemming from this quasi-equilibrium process at a de- 
lay of T^ d we determine an angular momentum transfer 
rate of <rf q d = 0.026^;^ /i B /ps considering M at 140 
K [30[| . In Tb the minimum of M is reached already af- 
ter 20 ps indicating a faster demagnetization, which is a 
consequence of the direct spin-lattice coupling (L = 3, 
c.f. Fig. [T]). The cooling mediated recovery of the initial 
magnetization is described by an exponential behavior 
during several 100 ps. In Gd diffusive cooling and slow 
demagnetization occur on similar time scales and lead 
to a plateau; in Tb cooling occurs after demagnetization 
and a recovery of M is observed at delays > 20 ps. 

Now three questions remain open, (i) What is the fast 
demagnetization time scale in Gd? (ii) Does Tb also 
show two distinct demagnetization time scales and if yes 
(hi) do both differ with respect to Gd? To answer these 
questions we employed fs x-ray pulses which we obtain by 
femtosecond slicing of the electron bunches in the storage 
ring 21, 3l[. Fig. H] confirms a clear reduction of M for 
both elements. In Gd we find after 3 ps a normalized 
XMCD signal of 0.7, identical to the level at which the 
slower demagnetization process sets in (inset in Fig. [3]). 
Employing the fs x-ray pulses we resolve the initial, fast 
demagnetization process in Gd. Also for Tb we find a 
sizeable drop of M within 2 ps (dashed areas in Fig. H]). 

To determine the characteristic time scales, the ps and 
fs time-resolved data have been fitted simultaneously by 
biexponential functions taking into account the different 
x-ray pulse durations (solid lines in Figs. [3]and|4]). For 
Tb we obtain r^ b = 8 ± 3 ps, which translates to an an- 
gular momentum transfer rate of aj^ = 0.29^q'o8 Mb/ps. 
We explain this process as being mediated by direct 
spin-lattice coupling under quasi-equilibrium conditions 
persisting at corresponding delays > 1 ps. This is 
much faster than cr^ d = 0.026 /xb/ps determined for 
the indirect interaction in Gd, which demonstrates that 
the direct spin-lattice coupling stemming from the non- 
spherical 4/ orbital distribution accelerates the demag- 
netization process in Tb. 

Next, we focus on the ultrafast demagnetization pro- 
cess. From our fits we determine within error bars identi- 
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FIG. 4: (color online) Time-dependent XMCD signals for Gd 
(top) and Tb (bottom) measured with fs x-ray pulses. Note 
the different time intervals. Solid lines depict biexponential 
fits determined by simultaneously fitting the fs and ps time- 
resolved data (cf. solid lines in Fig. [3J. A single exponential 
(dash-dotted line) for Tb yields an unsatisfactory fit. 



cal times rg d = 0.76 ± 0.25 ps and r e ^ b = 0.74 ± 0.25 ps. 
These times are shorter than reported for Gd/Fe mul- 
tilayers [l3| and similar to reports on TbFe alloys 
Since they are clearly longer than the pulse durations 
we rule out coherent processes promoted in Ref. Q. 
Note that our observations are likewise not compatible 
with demagnetization via superdiffusive spin transport 
[4|. The ultrafast, component of the demagnetization is 
50 % of the total loss in M and thus much too large 
to be explained by transport of the 5d valence electrons. 
We can furthermore exclude a mere transfer of the mag- 
netic moment from 4/ to 5d electrons because (i) the 
transferred moment is considerably larger than the va- 
lence electron spin polarization in both lanthanides. (ii) 
The 4/ — 5d states are coupled by intra-atomic exchange 
[32( ] . Depending on the coupling strength a transfer 
between 5d and 4/ magnetic moments or a concomi- 
tant change seems reasonable. We have performed time- 
resolved magneto-optical experiments, which probe pri- 
marily the valence band spin polarization. We find a 
reduction in the transient 5d magnetic signal concomi- 
tant with the XMCD signal measuring the 4/ magnetic 
moment. Considering that the optically excited elec- 
trons in Gd equilibrate with the crystal lattice during 
1 ps [26(| the ultrafast demagnetization occurs as long as 
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the system remains electronically excited. Our results 
demonstrate an efficient hot electron mediated momen- 
tum transfer to the lattice and determine the correspond- 
ing transfer rates of magnetic moment to the lattice to 
aj x b = 3.1±o | /i B /ps and c^ x d = l-5±g;][ /i B /ps. Com- 
pared to the quasi-equilibrium processes discussed above 
these rates are 10 and 50 times larger for Tb and Gd, re- 
spectively. We propose that the intraatomic 4/ — 5d ex- 
change interaction, about 100 meV [32], mediates this ac- 
celeration in the electronically excited state, as by means 
of the 4/ — 5d coupling a spin-flip scattering process in 
the conduction band affects the 4/ electrons as well and 
thereby drives the ultrafast demagnetization via indirect 
spin-lattice coupling [H[. 

Interestingly, in both 4/ elements the ultrafast process 
lasts longer than in the 3d transition metal ferromag- 
nets [H|, [111, [33| , although the momentum transfer rates 
are comparable, e.g. = 2.7 //b/ps for Ni [21|. This 
originates from the fact that the dominant part of M is 
generated by the 4/ electrons and is considerably larger 
for Gd and Tb than for Fe, Co, and Ni. As the magnetic 
moment carried by the valence electrons is considerably 
smaller than the 4/ one, several spin flips in the con- 
duction band of the lanthanide elements are required to 
obtain the same relative demagnetization as in a 3d ferro- 
magnet, where the magnetic moment resides completely 
in the conduction band and is much smaller. 

Finally, we compare our results with the model cal- 
culations of Ref. [5[ . In agreement we find laser-induced 
demagnetization of lanthanides on two time scales. How- 
ever, the orbital momentum of the 4/ shell cannot be ne- 
glected because we observe aj^ = 11 • cr^ d . Ref. [5] does 
not consider direct spin-lattice coupling, shown here to be 
essential, and predicts a figure of merit for the demagneti- 
zation time that is proportional to the ratio of Curie tem- 
perature and magnetic moment Tq / /i a t • Applying this to 
Gd (T c = 293 K, /i at = 7.55 /i B ) and Tb (T c = 225 K, 
/i a t = 9.34 jub) suggests that demagnetization in Gd is 
faster than in Tb by a factor of 1.6. Even within our 
conservative error bars our results cannot support this 
estimation; the faster demagnetization times coincide for 
Tb and Gd. Furthermore these times compare well with 
the time scale of electron-phonon equilibration [26]. We 
consider therefore that the time interval during which the 
faster demagnetization process is active is determined by 
electron-phonon interaction. A similar r ex for Gd and Tb 
is very plausible since their valence electron and crystal 
lattice are widely comparable. 

In conclusion magnetization dynamics of the 4/ mo- 
ments in Gd and Tb occurs on two timescales. The slow 
picosecond timescale is determined by the equilibrium 
spin-lattice coupling following the 4/ occupation. The 
fast femtosecond timescale is comparable for Gd and Tb 
and shows a pronounced enhancement of the valence- 
electron mediated indirect spin-lattice coupling. We ex- 



pect this mechanism to be operative also in the 3d fer- 
romagnets but hard to unravel due to the delocalized 
character of the magnetic moment. 
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